Protostellar Collapse with Various Metallicities 



Kazuyuki Omukai 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

ABSTRACT 

The thermal and chemical evolution of gravitationally collapsing protostellar clouds is 
investigated, focusing attention on their dependence on metallicity. Calculations are carried 
out for a range of metallicities spanning the local interstellar value to zero. During the time 
when clouds are transparent to continuous radiation, the temperatures are higher for those with 
lower metallicity, reflecting lower radiative ability. However, once the clouds become opaque, in 
the course of the adiabatic contraction of the transient cores, their evolutionary trajectories in 
the density-temperature plane converge to a unique curve that is determined by only physical 
constants. The trajectories coincide with each other thereafter. Consequently, the size of the 
stellar core at the formation is the same regardless of the gas composition of the parent cloud. 



Subject headings: ISM: clouds — ISM: molecules — molecular processes — stars: formation 
stars: Population II 
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1. Introduction 

In the standard theory of cosmic structure formation, galaxies are believed to have originated from 
primeval density fluctuations. Many authors have been studying the galaxy formation process in this 
context and have made remarkable progress in understanding such processes as the gravitational growth 
of the density fluctuation, its decoupling from the cosmic expansion and virialization, and the subsequent 
radiative cooling of the baryonic gas (e.g., Padmanabhan 1993). However, processes after the cooling of 
baryonic gas, namely, the transformation of a pregalactic cloud into a cluster of stars, are relatively poorly 
known. It is inevitably necessary to investigate them in studying the galaxy formation process since the 
very existence of stellar components is one distinct feature of galaxies. These processes can be viewed 
as successive fragmentations of a cloud into fragments, or in other words, protostellar clouds, and their 
contraction into stars (e.g., Hayashi 1984). 

Whether a cloud fragments or collapses into a single object depends on the the central flatness, that 
is, the axial ratio of the isodensity contour in the central region (Tsuribe & Inutsuka 1999). In particular, 
warm clouds do not fragment before adiabatic core formation, whereas clumps sufficiently more massive 
than the Jeans mass first collapse disklike and next fragment most likely into filamentary clouds (Miyama, 
Narita, & Hayashi 1987a, 1987b). The filamentary clouds fragment again into protostellar cloud cores after 
further contraction. Uehara et al. (1997) investigated the gravitational collapse of metal-free filamentary 
clouds using one-zone approximation and found that the minimum mass of fragments, or protostellar cores, 
is essentially Chandrasckhar mass, i.e., ~ IMq. Nakamura & Umemura (1999) confirmed Uehara et al. 
(1997)'s result by performing one-dimensional hydrodynamical calculations. In this paper, we will discuss 
the collapse of protostellar cores into stars. 

The distinction between the star formation in the galaxy formation epoch and that in present-day 
star-forming regions mainly resides in the difference of the composition of gas that stars are made from. 
Mctallicity, and therefore the gas to dust ratio, is lower for earlier star formation. In present-day star 
forming regions, protostellar clouds remain roughly 10K in a wide range of densities owing to efficient dust 
emission (e.g., Hayashi & Nakano 1965) in the course of contraction. 

On the other hand, many authors have studied the thermal and chemical evolution of collapsing 
primordial clouds with one-zone approximation (e.g., Matsuda, Sato, & Takeda 1969; Yoneyama 1972; 
Calberg 1981; Palla, Salpeter, & Stahler 1983; Izotov & Kolcsnik 1984; Lahav 1986; Puy & Signore 1997) 
or by hydrodynamical calculations (Matsuda et al. 1969; Villere & Bodenheimer 1987; Haiman, Thoul, & 
Loeb 1996; Omukai & Nishi 1998; Oliveira et al. 1998). The collapse of metal-free protostellar clouds, which 
is relevant for the first star formation in the universe, is induced by radiative cooling owing to molecular 
hydrogen lines, and their temperatures are about 1000K. 

Above are two extremes in composition of protostellar gas, namely, the former for the high-metallicity 
end and the latter for the low-metallicity end. 

Then, how does the collapse of slightly metal-polluted clouds proceed? Even in the galaxy formation 
epoch, significant metal enrichment can occur since the lifetime of massive stars (~ 10 6 yr) is shorter than 
the free-fall time of pregalactic clouds (<~ 10 8 yr). In fact, both Lya forest clouds and Population II stars 
are slightly polluted by heavy elements. Therefore, it is important to study star formation from slightly 
metal-polluted gas. In this paper, we aim to fill the gap between the protostellar collapse of primordial and 
present-day interstellar gas. 

This topic has been investigated in the pioneering works of Low & Lynden-Bell (1976) and Yoshii & 
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Sabano (1980). However, the following points remain to be fulfilled; in their works, (i) chemical reactions 
had not been solved, (ii) molecular coolants, which become important at high densities, had not been 
considered, and (iii) the evolution of clouds as dense as the stellar density (> 10 22 cm -3 ) had not been 
studied. With these points in mind, we investigate in this paper the thermal and chemical evolution of 
collapsing spherical clouds as a function of metallicity. 

The outline of this paper is as follows. In §2, we describe the method of our calculations. In §3, results 
of our calculations are presented. We summarize our work in §4. 



2. Model 
2.1. Basic Equations 

We consider a spherical cloud with mean metallicity Z . The helium concentration is assumed to be 
Vac = 0.0972 for all clouds.Q We assume that the dust-to-gas ratio is proportional to the mean metallicity, 
or in other words, that a fixed fraction of the heavy elements in the interstellar medium condenses into 
dust grains. We adopt the Pollack et al. (1994) model of grains in molecular clouds (see §2.2.2). For the 
case of local interstellar clouds (Z ~ IZq), the mass fraction of grain we used is 0.934 x 10~ 2 in the lowest 
temperature regime, and the gas-phase elemental abundances are yc — 0.927 x 10~ 4 , and t/q = 3.568 x 10 -4 , 
respectively, which correspond to 46% of oxygen and 72% of carbon depleted onto grains. These values are 
reduced proportionally for lower metallicity cases. We normalize the mean metallicity (and therefore also 
the dust-to-gas ratio and gas-phase metallicity in our model) relative to the local interstellar values Z\ oca \ 
and denote relative metallicity by z = Z/Z\ oca \. 

We calculate the time evolution of the central density, temperature, and chemical composition of the 
collapsing core. 

First, we assume that the dynamics is described by the free-fall relation, 

I - b < 3) 

where p is the density in the central region and the free-fall time is 



'■■ 3 ^m,- (4 » 



The thermal evolution is followed by solving the energy equation 

de d , 1 

di ~ ~ P di^p 



''' vi{-)-^\ (5) 



The concentration of element X is defined by 

yx = nx/n, (1) 
where n and nj are the number densities of hydrogen nuclei and nuclei of element X. Similarly, we write for each atomic, 
molecular, or ionic species 

y(x)=n(x)/n, (2) 
where n(x) is the number density of species x. Note {/(H2) = 1/2 for fully molecular gas. 
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where 

e= (6 

7ad - 1 fim H 

is the specific internal energy, 

(7) 

is the pressure for an ideal gas, 7 a( j is the adiabatic exponent, T is the temperature, fj, is the mean molecular 
weight, toh is the mass of hydrogen nucleus, and £( not ) is the net energy loss rate per unit mass (see §2.2). 

We neglect any effect owing to rotation or magnetic fields for simplicity. In this case, the actual collapse 
is expected to proceed like the Penston-Larson similarity solution (Penston 1969; Larson 1969). According 
to this solution, the cloud consists of two parts, that is, the central core region, which has flat density 
distribution, and the envelope, where the density decreases outward as oc r~ 2 . The size of the central flat 
region is roughly given by the local Jeans length Aj = Trc s /y/Gp in the core. Since we focus on the evolution 
of the central region, optical depth is estimated by that across one local Jeans length: 

t v = n v p\ } . (8) 

This formulation of Jeans length shielding is the same procedure as that of Low & Lynden-Bell (1976) and 
Silk (1977), although they introduced it as a result of successive fragmentation into Jeans mass clouds, not 
because of the Penston-Larson collapse. 



2.2. Cooling/Heating processes 

In addition to compressional heating, we include cooling owing to (i) atomic and molecular line 
radiation, (ii) energy transfer between the gas and the dust grains, which will be emitted as infrared 
radiation from the grains, (iii) continuous radiation from the gas, and (iv) cooling and heating associated 
with chemical reactions. Then 

£ (nct) = £] 

1110 cont ~f~ ^chemj (9) 

where the terms in the right-hand side of the equation correspond to the cooling rates owing to the processes 
(i)-(iv) above in the same order. Although we treat clouds with metallicity up to ~ IZq, we are concerned 
mainly with the lower metallicity environments that are relevant to the epoch of galaxy formation. In this 
paper, we do not try to reproduce the present Galactic environment and neglect any external radiation (e.g., 
UV radiation, cosmic rays, cosmic background radiation, etc.) for simplicity. Then, the considered heating 
sources are the compressional work and H2 formation. As will be mentioned below, lower metallicity clouds 
tend to be warmer. Therefore, the compressional heating, which is proportional to the temperature of the 
cloud, becomes more important if the external heating rate is constant. 



2.2.1. Line Cooling 

As line cooling agents, we include atomic fine-structure transitions of CI, CII, and 01 and molecular 
rovibrational transitions of H 2 , CO, OH, and H 2 0: 



Aine = £ci + £cil + + £h 2 + £cO + ^OH + £h 2 



(10) 
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The line cooling rate of species x is given by 

C x = - n{x,i)Aijajhuij, (11) 
P (i-j) 

where n(x, i) is the population density of species x in level i, Aij is the spontaneous transition probability, 
tij is the escape probability and hvij is the energy difference between levels % and j. The population density 
n(x,i) is obtained from a solution of the equation of statistical equilibrium 

n n 

n{x,i)^Rij =^2n(x,j)Rji, (12) 

where n is the total number of lines included and 

„ = f Aijtij + dj for i > j 
13 I Cij for i < j, 

ignoring the external radiation. Here Cij is the collisional rate from level i to level j. Following Takahashi, 
Hollenbach, & Silk (1983), we use the escape probability for the case that the velocity is proportional to the 
radius in the central region: 

' 1 - c~ Ti j \ 

(14) 



lij 

The optical depth averaged over the line and the continuum optical depth are given by 

A- c 3 

T *J = i n ( x J)9i/9j ~ n(x,i)]l sh /(2Av D ), (15) 

and 

Tcont = («gr + K gas )p\j, (16) 

respectively, where ijt is the multiplicity factor, gt is the statistical weight of level i, A«d is the velocity 
dispersion, and l sh is the shielding length (Takahashi et al. 1983): 

dv 

l sh = min(As t h, Aj); As th = 2Aw D /( — ) = 6Aw D t ff , (17) 

dr 

where we have used the relation for the homogeneous collapse v = r/3ts and, in the case of the small 
velocity gradient, have assumed Jeans length shielding as introduced in §2. As a source of the velocity 
dispersion, we consider only the thermal motion of atoms and neglect microturbulent motions. Then 



Av„ = J^, (18) 

where fi x is the molecular weight of species x. The continuum Planck mean opacity of dust K gr and that of 
gas /£gas are adopted from Pollack et al.(1994) and Lenzuni, Chernoff, & Salpeter (1991), respectively (see 
§2.2.2 and 2.2.3). 

In general, the level population and the escape probabilities depend on each other. We then need 
an iterative procedure to find their consistent solution except for the case of H 2 . For H 2 lines, we need 
not iterate to find the consistent set of the level populations and the escape probabilities, since in our 
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calculations, clouds become opaque to H2 lines only at the density exceeding the critical density, and in 
such a case, molecules populate following a Boltzmann distribution regardless of the escape probabilities. 

The related parameters of transitions are given in Hollenbach & McKee (1989) for CI, CII, and 01, in 
Hollenbach & McKee (1979;1989) for OH and H 2 Q and in McKee et al. (1982) for CO. We computed 
the population of H2 following the procedure of Hollenbach & McKee (1979) using the collision coefficient 
given in Hollenbach & McKee (1989). For H2, we considered the first three vibrational states (v = — 2) 
with rotational levels up to J = 20 in each vibrational state. We take into account only the ground 
vibrational level for other molecular lines since cooling rates owing to vibrational transitions are nearly 
always dominated by grain cooling or H 2 cooling (Hollenbach & McKee 1979). 



2.2.2. Gas- Grain Heat Transfer 

We adopt the Pollack et al. (1994) model of grains in molecular clouds. In their model, silicates, 
organics, troilite, metallic iron, and water ice constitute the most abundant grain species and the dominant 
sources of opacity in molecular cloud cores. The Planck mean opacity due to the grains in the case of mean 
metallicity Z = Z\ oca \ is 

Kgr = 4.0 x 10~ 4 T 2 r cm 2 g-\ (19) 

in T gr <; 50K, where T gr is the effective grain temperature (see below). The mass fraction of grain is 
0.934 x 10~ 2 below the water vaporization temperature (100-200 K, depending on the density), and that 
of each component; water ice, volatile organics, refractory organics, troilite, orthopyroxene, olivine and 
metallic iron are 1.19 x 10~ 3 , 6.02 x 10~ 4 , 3.53 x 10~ 3 , 5.69 x 10~ 4 , 7.33 x 10~ 4 , 2.51 x 10~ 3 , and 2.53 x 10~ 4 , 
respectively (from Table 2 of Pollack et al. 1994). For the grain size distribution 



1 a (0.005um < a < lam) 

" (a)0C ^ (l M m< a <5 M m) (20) 



(Pollack et al. 1994) that is based on that derived by Mathis, Rumpl, & Nordsieck (1977), the energy 
transfer rate from the gas to the dust grains per unit mass is given by 

£ gr = 1.1 x 10 5 n(^)(— ^) 1/2 [1 - 0.8exp(-75K/T)](T- T gr ) (21) 

(Hollenbach & McKee 1979), where T gr is an effective grain temperature, which is determined by the energy 
balance for dust grains; 

£ gr = 4aT 4 I .K gI ./3 collt , (22) 

where n gI is the Planck mean opacity owing to the dust grains and / gr /p gr is the total volume of dust 
grains per unit mass of gas and / gr /p gr = 5.3 x 10 -3 at the lowest temperatures (i.e., below the water ice 
vaporization temperature). These quantities are taken from tables in Pollack et al.(1994). 

The continuum energy transport rate /3 C ont decreases as r~ ^ t owing to radiative diffusion in the 
optically thick case (e.g., Masunaga et al. 1998). Then, 

Axmt = min(l,T~ 2 lt )- (23) 



2 The OH and H2O collision cross sections and H2O multiplicity factor are given in Hollenbach & McKee (1989). We adopted 
other parameters from Hollenbach & McKee (1979). 
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2.2.3. Continuum of Gas 

Using the Planck mean opacity of gas K gas , the cooling rate owing to the gas continuum is given by 

£cont = 4oT 4 Kg as /3 cont . (24) 

We take the continuum Planck mean opacity for metal-frcc gas from Lcnzuni ct al. (1991), which includes all 
the important continuum processes, specifically, bound-free absorption by H° and H~ ; free-free absorption 
by H°, H~, H 2 , H^~, Hj, Hg , He , He - ; photodissociation of H 2 , and H^ by thermal radiation; Raylcigh 
scattering by H°, H 2 , He ; Thomson scattering by e~; and collision- induced absorption by H 2 due to 
collisions with H 2 , He , and H . Among the above processes, the most important cooling mechanism is H 2 
collision-induced continuum, which dominates the cooling at about <~ 10 16 cm~ 3 for low metallicity (i.e., 
< 1CT 6 %) clouds. 



2.2.4- Chemical Cooling/Heating 

Following Hollenbach & McKee (1979), we assume the heat deposited per a formed molecular 
hydrogen as 0.2 + 4.2(1 + n cr /n) _1 cV for H 2 formation on grain surfaces (reaction 23 in Appendix), 
3.53(1 + n cr /n) _1 eV for H 2 formation by H~ process (reaction 8), 1.83(1 + n cr /n) _1 cV for H 2 formation 
by H^ process (reaction 10), and 4.48(1 + n cr /n) _1 eV for H 2 formation by the three-body reactions 
(reactions 19 and 20), where 

in 6 T-V 2 

cm" 3 . (25) 



1.6y(H)cxp[-(400/T) 2 ] + 1.4j/(H 2 )exp[-12000/(T + 1200)] 

Collisional dissociation and ionization absorb the same energy as the binding energy, that is 4.48 eV per H 2 
dissociation, 13.6 eV per H ionization, 24.6 eV per He ionization, and 54.4 eV per He+ ionization. 



2.3. Chemical reactions 

We solve nonequilibrium chemistry involving the four elements H, He, C and O, that contains the 
following 45 species: H, H 2 , e", H+, H+, H+, H", He, He+, He++, HcH+, C, C 2 , CH, CH 2 , CH 3 , CH 4 , C+, 
C+, CH+, CH+, CH+, CH+, CH+ O, 2 , OH, CO, H 2 0, HCO, 2 H, C0 2 , H 2 CO, H 2 2 , 0+, 0+, OH+, 
CO+, H 2 0+, HCO+, 2 H+, H 3 0+, H 2 CO+, HCOj and H 3 CO+. Chemical reactions included arc listed in 
Appendix. The reactions of H and He chemistry are mainly taken from Abel et al. (1997) and Galli & Palla 
(1998). Other important H 2 forming processes included are the three-body reactions (Palla et al. 1983) and 
the reaction on the surfaces of the dust grains (Tielens & Hollenbach 1985). In addition, we supplement H, 
He, C and O chemical reactions involving the above species from Millar, Farquhar, & Willacy (1997). 



3. Results 

In this section, we present the results obtained by the method described in §2 and discuss our analysis. 

Figure 1 displays the evolutionary trajectories of protostellar clouds whose metallicities are 
z = 0, 10~ 6 , 10~ 4 , 10~ 2 , 1. Here we set the initial condition to be T = 100K,n = lc m - 3 ,y(e) = 
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1x10 4 , y(H 2 ) = 1 x 10 6 . All the carbon is assumed to be in the form of CII, while oxygen is 01 at the 
beginning. 

The evolutionary trajectories for two other initial conditions are shown in Figure 2, as well as the same 
curves in Figure 1. We can see from Figure 2 that trajectories of clouds with a fixed composition converge 
rapidly toward the dense region for any initial conditions (e.g., Hayashi & Nakano 1965; Low & Lynden-Bell 
1976). We have also tested the sensitivity to initial chemical compostions for two cases - (i) where the 
ionization degree y(e) = 1 x 10~ 3 and (ii) all the carbon is assumed to be CI,- and have found the results 
similar to those shown in Figure 2. Hereafter, we discuss only on the case presented in Figure 1. 

In general, as seen in Figure 1, the temperatures of lower metallicity clouds are higher because of their 
lower radiative cooling ability (i.e., lower radiative cooling rate for the same temperature and density) as 
long as the clouds are transparent to continuum and continue to collapse dynamically, owing to the efficient 
radiative cooling, in other words, during what is known as the "first collapse" stage. 

On account of the total lack of metals and grains, the only cooling agent in the temperature range 
below <~ 10 4 K for primordial (z = 0) clouds is rovibrational transitions of molecular hydrogen (e.g., Matsuda 
et al. 1969). In the primordial gas, H 2 is formed mainly by the H~ process, 

H + er -» H-+ 7 ; (26) 
H + H- -» H 2 +e-, (27) 

(Saslaw & Zipoy 1967) until the density reaches ~ 10 8 cm -3 , where the three-body reactions 

3H -» H 2 + H (28) 

2H + H 2 -> 2H 2 (29) 

become efficient (Palla et al. 1983). H 2 line emission via electric quadrupole transitions is the only 
efficient cooling agent until the number density exceeds ~ 10 14 cm -3 , where H 2 continuum emission via 
collision-induced dipole transitions begins to overwhelm. The drop in temperature at n ~ 10 16 cm~ 3 is 
due to the quadrature dependence of the collision-induced emission coefficient on the density. The cloud 
becomes optically thick to the continuum at 10 16 cm -3 . (See Fig. 3 a) 

The gravitational contraction of metal-free protostellar clouds has been investigated with Omukai 
& Nishi (1998) by detailed hydrodynamical calculations assuming spherical symmetry. Our evolutionary 
trajectory for the primordial cloud agrees well with those of Omukai & Nishi (1998) at the number density 
n J> 10 10 cm~ 3 . At lower densities, the temperature of ours is higher owing to the influence of Omukai & 
Nishi's (1998) initial condition. They started the calculation from a cloud in hydrostatic equilibrium at 
n ~ 10 6 cm~ 3 . As long as the density is not so much higher than the initial value, the collapse is slower 
than the free-fall rate, which we have assumed in this paper. As a result, the temperature is lower because 
of the lower compressional heating. 

The evolution of the cloud with z = 10 -6 is essentially the same as that of the primordial cloud, 
although there is minor deviation caused by extra cooling owing to water molecules and dust grains. (See 
Fig. 3 b) The dust thermal emission dominates the cooling just before the H 2 continuum becomes effective. 
However, it is only temporary and the grains rapidly sublimate at n~ 10 14 cm~ 3 . The most refractory grain 
composition that contributes to the Planck mean opacity is iron, which sublimate at about 1200K at this 
density (Table 3 of Pollack et al. 1994). The slight amount of dust grains in the cloud with z = 10~ 6 does 
not contribute to H 2 formation significantly, as seen in Figure 4. 
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For the cloud with z = 10~ 4 , H2, which is formed mainly on the grain surfaces instead of the H~ 
process, is the major coolant at low densities (< 10 4 cm~ 3 ). Water molecules, which are mainly formed at 
~ 10 6 cm -3 by the reaction 

H 2 + OH^H 2 + H, (30) 

dominate the cooling in the range 10 5 cm~ 3 < n < 10 10 cm~ 3 . In present-day molecular cloud cores, usually 
neither H2O nor OH plays a key role in the thermal balance since the temperature there is lower than the 
excitation energy of these molecular species (Neufeld, Lepp, & Melnick 1995). However, in lower metal 
environments, the temperature is higher because of the lower radiative cooling. In that case, H 2 can be 
a principal cooling agent. In the cloud with z = 10~ 4 , grain surface reactions exceed the H~ process as a 
mode of H 2 formation, although hydrogen does not become fully molecular until the three-body reaction 
begins to work. At 10 10 cm -3 , the grain thermal emission dominates the cooling and causes the temperature 
to drop to about 100K. The cloud becomes optically thick to the grain emission at about 10 13 cm~ 3 . (See 
Fig. 3 c) 

For the cloud with z > 10~ 2 , cooling is dominated by atomic lines and CO lines in lower densities (that 
is, n < 10 6 cm -3 for z = 1CP 2 , and n < 10 3 cm~ 3 for z = 1, respectively), and by the dust thermal emission 
in higher densities. (See Fig. 3 d, e) The temperature of our z = 1 cloud drops far less than 10K, which is 
appropriate in present-day protostellar clouds. This is merely because we have ignored external radiation, 
and the results should not be taken seriously as a model of present-day molecular cloud cores. 

Clouds with metallicity z > 10~ 4 become opaque to the thermal radiation of dust grains before the 
grain vaporization occurs. Once the cloud becomes opaque, it begins to contract adiabatically since the 
radiative cooling rate drops rapidly. According to some hydrodynamical simulations (e.g, Larson 1969; 
Masunaga et al. 1998), a transient core in hydrostatic equilibrium (called a "first core") forms after the 
cloud has more contractions. The hydrostatic core continues the adiabatic contraction by accreting the 
envelope matter until the temperature reaches ~ 2000 K, where molecular hydrogen begins to dissociate. 

Clouds with lower metallicity (i.e., z < 10~ 6 ) become opaque to the H 2 continuum instead of the dust 
thermal emission, as described above. Since the temperature at that time is already near the dissociation 
value, molecular hydrogen begins to dissociate just after the cloud becomes opaque. As a result, no transient 
hydrostatic core forms in these cases as has been demonstrated by the hydrodynamical calculations of 
Omukai & Nishi (1998). 

When the ratio of specific heats T = dlogp/dlogp exceeds the critical value 4/3, corresponding to a 
gradient of 1/3 in the density-temperature (n — T) plane (i.e., Figure 1), the bulk motion of the cloud is 
decelerated and the free-fall assumption is no longer valid. However, the trajectories in the n — T plane 
are not altered since the clouds contract almost adiabatically anyway in the opaque stage because of the 
rapidly declining radiative cooling rate (e.g., Larson 1969; Narita, Nakano, & Hayashi 1970). 

In the course of the adiabatic contraction of the transient cores, all the evolutionary trajectories 
converge to a certain line (Fig. 1, dashed line) in spite of their different composition and histories. We 
discuss the reason for the convergence in the following. The trajectories in the n — T plane just prior to 
when the clouds become opaque are determined by thermal balance between the compressional heating and 
the radiative cooling. The compressional heating rate (per unit mass) is given by 

Qc °^ = - p Jt { )) = |' (31) 

where p is the pressure and c s is the sound speed in the central region of the cloud. The radiative cooling 
at that time is dominated by continuum for all clouds, namely, the dust thermal emission for clouds with 
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relative metallicity z > 10 4 , or the H2 continuum emission for those with lower metallicity. The radiative 
cooling rate is given by 

£rad = ^GT\ (32) 

where k b is the Planck mean continuum opacity. In the equation above, we have used the relation that 
the grain temperature is the same as the gas temperature, T gr = T, which is valid in a cloud so dense that 
it becomes opaque to the continuum. Note that expression (p2|) holds for any clouds that cool mainly 
by continuum by choosing an appropriate value of opacity, where the opacity is an increasing function of 



relative metallicity z. Equating (31) and (|32|), for each value of k b (z), we obtain a trajectory just before 



the cloud becomes opaque (we call it Tr(« B )); 



Tv(k b ) : f = An B aT 4 . (33) 



As the cloud of opacity k,b collapses along the trajectory represented by (|33|), the central part of the cloud 
becomes opaque to the continuum at the point 

P(k b ) : tj(k b , n, T) = K B p7r— ^= = 1, (34) 
where we have used the assumption that the size of the central region is approximately Aj. Eliminating k b 



in equations ( |33| ) and (34), we obtain the locus L of points P for variable k b ; 

L: r = <s^> ,/5 " 2/5 - < 35 » 

On the other hand, the trajectory of the cloud of opacity k b in the opaque stage, 0(k b ), is the line that 
starts from the point P(k b ) and whose gradient is approximately equal to the value for adiabatic contraction 
7 a d — 1. The adiabatic exponent 7 a d for molecular hydrogen is 7/5 in the temperature range T ^ 200K, 
while it is 2/3 in the lowest temperature. The gradient of the trajectory 0(k b ) in the range T ;> 200K is 
2/5. This value happens to be the same as the gradient of the line L. Thus, the trajectories in the opaque 
stage O all coincide with the line L, which can be written only by physical constants. When the clouds 
become opaque (i.e., at the point P), their specific entropy has the same value for all clouds regardless of 
their metallicity. 

As a cloud climbs up along line L, the temperature reaches the dissociation value (~ 2000K) at 
n~3x 10 16 cm~ 3 . The ratio of specific heats are reduced below the critical value 4/3. Then the cloud 
begins to collapse dynamically again, or in other words, it begins the second collapse (e.g., Larson 1969). 
The trajectories in the second collapse stage are also unique because the specific entropy has the same value 
for all clouds regardless of their metallicity. When most hydrogen molecules have been dissociated, the ratio 
of specific heats rises above the critical value 4/3 again. After the central part of the cloud contracts almost 
adiabatically to some extent, a hydrostatic core forms. This core is referred to as the stellar or second 
core in literature. We cannot accurately estimate the size of the stellar core from our simple one-zone 
treatment. However, hydrodynamical calculations tell us that the number density and mass of the stellar 
core at its formation time are on the order of 10 22 cm -3 and 10 -3 Mq respectively, both for present-day 
(Larson 1969) and primordial protostars (Omukai & Nishi 1998). This value is expected to be universal, 
namely, independent of metallicity since the equations of state, or trajectories in the n — T plane, are the 
same in the second collapse stage. 



The stellar core, or protostar, grows in mass by accretion of the envelope material. The temperature 
rises accordingly and eventually the core will become an ordinary star. Consequently, the final mass of stars 
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is determined by the subsequent accretion onto the core, although the size of the protostar at the time of 
formation is the same. Stahler, Shu, & Taam (1980) argued that a rough estimate of the protostellar mass 
accretion rate M can be obtained from the relation 

M ~ 4/G, (36) 

where c s is the isothermal sound speed in the initial protostellar cloud. Therefore, the mass accretion rate 
is higher for a lower metallicity cloud because of the higher temperature of the protostellar cloud (Stahler, 
Palla, & Salpeter 1986). The higher accretion rate and lower opacity (then smaller radiation force) may 
result in the higher mass of formed stars for a lower metallicity cloud (e.g., Wolfire & Cassinelli 1987). 
However, to address these issues fully, further investigations are needed. 

4. Summary 

We have investigated the thermal and chemical evolution of collapsing protostellar clouds with different 
metallicities. The varied range of metallicity spans the local interstellar value (~ IZq) to zero. 

The evolution of the clouds is summarized as follows. While the clouds are transparent to continuous 
radiation, the temperature of clouds with lower metallicity is higher since their radiative cooling rates are 
lower for the same density and temperature. However, after the clouds become opaque and begin adiabatic 
contraction, their evolutionary trajectories converge to a line that is determined only by physical constants. 
Thereafter, the trajectories coincide with each other regardless of their metallicity Consequently, the 
physical dimension of the stellar core at the time of formation is the same for clouds with any composition. 

We have also discussed analytically the reasons for the convergence. 

The author wishes to express his cordial thanks to Humitaka Sato and Naoshi Sugiyama for their 
continual interest and advice, to Hiroshi Koyama, Ryoichi Nishi, and Toru Tsuribe for useful discussions, 
and to the referee for improving this manuscript. 
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Table 1. Chemical Reactions Included 



Number Reaction Rate Coefficient a Reference 

1 H + c -» H+ + 2c fci = cxp[-32. 71396786 + 13.5365561n(T(eV)) A97 

-5.73932875(ln(T(eV))) 2 
+1.56315498(ln(T(eV))) 3 
-0.2877056(ln(T(cV))) 4 
+3.48255977 x 10- 2 (ln(T(eV))) 5 
-2.63197617 x 10- 3 (ln(T(cV))) 6 
+1.11954395 x 10- 4 (ln(T(cV))) 7 
-2.03914985 x 10- 6 (ln(T(cV))) 8 ] 

2 H++c^H + 7 k 2 = cxp[-28. 6130338 - 0.72411256(ln(T(eV))) A97 

-2.02604473 x 10- 2 (ln(T(eV))) 2 
-2.38086188 x 10- 3 (ln(T(eV))) 3 
-3.21260521 x 10- 4 (ln(T(cV))) 4 
-1.42150291 x 10- 5 (ln(T(eV))) 5 
+4.98910892 x 10- 6 (ln(T(cV))) 6 
+5.75561414 x 10- 7 (ln(T(eV))) 7 
-1.85676704 x 10- 8 (ln(T(eV))) 8 
-3.07113524 x 10- 9 (ln(T(cV))) 9 ] 

3 He + c -» Hc+ + 2e fc 3 = cxp[-44.09864886 + 23.91596563(ln(T(eV))) A97 

-10.7532302(ln(T(eV))) 2 
+3.05803875(ln(T(eV))) 3 
-0.56851189(ln(T(eV))) 4 
+6.79539123 x 10- 2 (ln(T(cV))) 5 
-5.00905610 x 10- 3 (ln(T(eV))) 6 
+2.06723616 x 10- 4 (ln(T(cV))) 7 
-3.64916141 x 10- 6 (ln(T(eV))) 8 ] 

4 Hc++c^Hc + 7 fc 4 = 3.925 x 10- 13 T(cV)-°- 6353 A97 

+ 1.544 x 10- 9 T(eV)- 15 
[3.0 x 10- 1 exp(-48.596/T(eV)) +exp(-40.496/T(eV))] 

5 He+ + e -> Hc++ + 2c fc 5 = cxp[-68. 71040990 + 43.93347633(ln(T(cV))) A97 

-18.4806699(ln(T(eV))) 2 
+4.70162649(ln(T(eV))) 3 
-0.76924663(ln(T(eV))) 4 
+8.113042 x 10- 2 (ln(T(eV))) 5 
-5.32402063 x 10- 3 (ln(T(cV))) 6 
+1.97570531 x 10- 4 (ln(T(cV))) 7 
-3.16558106 x 10- 6 (ln(T(cV))) 8 ] 

6 Hc+++c^Hc++ 7 fc 6 = 3.36x lO-^T- 1 / 2 ^/^ 3 )- - 2 !^^/^ 6 ) - 7 ]- 1 A97 

7 H + e^H-+7 fc 7 = 1.4 x 10- 18 T°- 928 exp(-r/1.62 x 10 4 ) GP98 

8 H" + H -> H 2 + e fc 8 = 4.0 x 10- 9 T-° 17 (T > 300) GP98 

fc 8 = 1.5 x 10~ 9 (T < 300) 

9 H + H+ — > + 7 k 9 = iq[-W-M-i-5M1o Sio t+i.iis(1os 10 t) 2 -o.1269(1o Sio t) 3 ] GP98 



Tabic 1 — Continued 



Number 



Reaction 



Rate Coefficient a 



Reference b 



fO 
ff 

12 
13 



14 



22 
23 



24 
25 



H++H- 
H 2 +H+ 
H 2 +e^ 
H 2 + H- 



♦ H 2 + H+ 
-^H++H 
2H + e 
3H 



H" 



H + 2e 



2H — > H+ + e + H 
2H + grain -> H 2 



He+ + H 2 
H+ + He - 



-> H+ + H + He 
HeH+ + H 



15 


H" + H+ -» 2H 


fcl5 


16 


H- + H+ -> H+ + e 


&16 






fcl6 


17 


H+ + c -» 2H 


kn 


18 


H+ + H- -» H 2 + H 


kis 


19 


3H -» H 2 + H 




20 


2H + H 2 -» 2H 2 


k 2 o 


21 


2H 2 — > 2H + H 2 


k 2 i 



fcio = 6.4 x 10~ 10 

fen = 3.0 x 10- 10 cxp(-2.1050 x 10 4 /T) 
fci2 = 4.4 x 10- 10 T a35 exp(-1.02 x 10 5 /? 1 ) 

*i3 = °*£ 

fc L = 1.12 x 10- 10 cxp(-7.035 x 10 4 /T) 
fc H = 6.5 x 10- 7 T-V2 

cxp(-5.2 x 10 4 /T)[1 - cxp(-6000/T)] 
a = (1 + n/n cr ) _1 

log 10 (n cr ) = 4.0 - 0.4161og 10 (T/10 4 ) 
-0.327(log 10 (T/10 4 )) 2 
few = cxp[-18. 01849334+ 2.3608522(ln(T(eV))) 
-0.28274430(ln(T(cV))) 2 
+ 1.62331664 x 10- 2 (ln(T(eV))) 3 
-3.36501203 x 10- 2 (ln(T(eV))) 4 
+ 1.17832978 x 10~ 2 (>(T(cV))) 5 
-1.65619470 x 10- 3 (ln(T(eV))) 6 
+ 1.06827520 x 10- 4 (ln(T(cV))) 7 
-2.63128581 x 10~ 6 (>(T(eV))) 8 ] 
5.7 x lO- 6 ^ 1 / 2 

+ 6.3 x 10~ 8 - 9.2 x lO" 11 ! 11 / 2 + 4.4 x 10~ 13 T 
6.9 x io- 9 T-°- 35 (T < 8000) 
9.6 x 10- 7 T-°- 9 (T > 8000) 
2.0 x 10- 7 T-V 2 
2.30 x 10- 7 (T/300)- 1 / 2 
5.5 x lO" 29 ! 1 - 1 
kw/8 

H K L 

k L = 1.18 x 10- 10 cxp(-6.95 x 10 4 /T) 
fc H = 8.125 x IQ-^T- 1 / 2 

exp(-5.2 x 10 4 /T)[1 - cxp(-6.0 x 10 3 /T)] 
a = (1 + n/n^y 1 

log 10 (n cr ) = 4.845- 1.31og 10 (T/10 4 ) 
+1.62(log 10 (T/10 4 )) 2 
k 22 = 1-7 x 10- 4 /ci 
fc 23 = 6.0 x 10- 17 (T/300K) 1 / 2 / a 
x [1 + 4.0 x 10- 2 (T + T gr ) 1 / 2 

+2.0 x 10~ 3 T + 8.0 x lO^T 2 ]- 1 
x Z/Z\ oca \ 

/a = [1 + exp(7.5 x 10 2 (l/75 - 1/T^))]- 1 
k 24 = 3.70 x 10~ 14 exp(-35/T) 
k 2b = 1.30 x 10- 10 



GP98 
GP98 
GP98 

SK87 



A97 
GP98 

GP98 

GP98 

M97 

PSS83 

PSS83 

SK87 

PSS83 (Afe) 



PSS83 
TH85 



M97 
M97 
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Table 1 — Continued 



AT 1_ 

Number 


Reaction 








Rate Coefficient a 




Tl c b 

Rctcrence 


26 


H^ + H 2 - 


-> Rt +H 






2.08 X 


10" 






M97 


27 


H+ + H- - 


-» 2H 2 


IV Z I 




2.30 x 


10" 


- 7 (T/300)-°- 50 




M97 


28 


He+ + H - 


■+ H+ + He 


Too 

"'28 




1.90 x 


10- 


-15 




M97 


29 


Hc+ + H- 


-> H + Hc 






2.30 x 


10" 


- 7 (T/300)-°- 50 




M97 


30 


He+ + H 2 


-» H+ + He 


kin 




7.20 x 


10" 


-15 




M97 


31 


HeH+ + H 


-y H+ + He 


hn 




9.10 x 


10" 


-10 




M97 


32 


HcH+ + H 


2 -» H+ + He 


kio 




1.50 x 


10" 


-9 




M97 


33 


H+ + e^ 


H 2 +H 


kii 

nj oo 




1.50 x 


10" 


- 8 (T/300)-°- 50 




M97 


34 


H+ + e^ 


3H 


kiA 

n -o4 




1.50 x 


10" 


- 8 (T/300)-°- 50 




M97 


35 


HeH+ + e 


-» H + He 


Kr 

"'00 




2.00 x 


10" 


- 7 (T/300)-°- 50 




M97 


36 


H + CH^ 


C + H 2 


kia 
"'Jo 




4.98 x 


10" 


-11 




M97 


37 


H + CH^ 


C + 2H 


ki7 




6.00 x 


10" 


- 9 exp(-40200/T) 




M97 


38 


H + CH 2 - 


■+ CH + H 2 


kiQ 




2.70 x 


10" 


-10 




M97 


39 


H + CH 3 - 


-> CH 2 + H 2 


kin 




1.00 x 


10" 


- 10 exp(-7600/T) 




M97 


40 


H + CH 4 - 


-> H 2 + CH3 


ft! /in 




5.82 x 


10" 


13 (T/300) 3 00 exp(- 


-4045/T) 


M97 


41 


H + OH- 


•H 2 + 


A* /i 1 
"-41 




7.00 x 


10" 


- 14 (T/300) 2 - 80 exp(- 


-1950 /Tl 


M97 


42 


H + OH- 


• + 2H 


k,Ai 




6.00 x 


10" 


' 9 exp(-50900/T) 




M97 


43 


H + H 2 0- 


-+ OH + H 2 


k,Ai 




6.83 x 


10" 


- 12 (T/300) 1 - 60 exp(- 


-9720/T) 


M97 


44 


H + H 2 0- 


OH + 2H 


' v 44 




5.80 x 


10" 


- 9 exp(-52900/T) 




M97 


45 


H + C 2 -> 


CH + C 


r *-4o 




4.67 x 


10" 


- 10 (T/300)°- 50 exp(- 


-30450 /T) 


M97 


46 


H + COh 


C + OH 


,v 40 




1.10 x 


10" 


10 (T/300)°- 50 exp(- 


-77700/T) 


M97 


47 


H + H 2 CO 


-» HCO + H 2 


At /i ? 

,v 4 1 




1.50 x 


10" 


- n (T/300) 1 05 exp(- 


-1650/T) 


M97 


48 


H + 2 -y 


OH + O 


'"4o 




3.30 x 


10" 


- 10 exp(-8460/T) 




M97 


49 


H + 2 -» 


20 + H 






6.00 x 


10" 


" 9 exp(-52300/T) 




M97 


50 


H + 2 H - 


-» H 2 + O 


"-50 




5.00 x 


10" 


n exp(-866/T) 




M97 


51 


H + 2 H- 


-+ H 2 + 02 






7.10 x 


10" 


- n exp(-710/T) 




M97 


52 


H + 2 H- 


-> 20H 






2.80 x 


10" 


10 exp(-440/T) 




M97 


53 


H + H 2 2 


-» 2 H + H 2 


k^i 




2.80 x 


10" 


- 12 exp(-1890/T) 




M97 


54 


H + C0 2 - 


^ OH + CO 


k^A 




2.50 x 


10" 


' 10 exp(-13300/T) 




M97 


55 


C + H 2 -» 


CH + H 


k^ 




6.64 x 


10" 


- 10 exp(-11700/T) 




M97 


56 


C + OH^ 


CH + O 


k^a 




2.25 x 


10" 


- n (T/300) - 50 exp(- 


-14800/T) 


M97 


57 


C + CO^ 


c 2 + o 


&57 




2.94 x 


10" 


11 (T/300)°- 50 cxp(- 


-58025/T) 


M97 


58 


+ H 2 -» 


OH + H 


&58 




3.43 x 


10" 


- 13 (T/300) 2 67 exp(- 


-3160/T) 


M97 


59 


O + CH^ 


• OH + C 


&59 




2.52 x 


10" 


- n exp(-2381/T) 




M97 


60 


+ CH 2 - 


^ CH + OH 


^60 




4.98 x 


10" 


- 10 exp(-6000/T) 




M97 


61 


+ CH 2 - 


-» HCO + H 


fcei 




1.44 x 


10" 


- n (T/300) - 50 exp(- 


-2000/T) 


M97 


62 


+ CH 4 - 


-» CH 3 + OH 






8.32 x 


10" 


- 12 (T/300) 156 exp(- 


-4270/T) 


M97 


63 


+ H 2 - 


-» 20H 


&63 




1.85 x 


10" 


- n (T/300) - 95 exp(- 


-8571/T) 


M97 


64 


+ H 2 CO 


- -» HCO + OH 


fc 64 




1.78 x 


10" 


- n (T/300) - B7 exp(- 


-1390/T) 


M97 


65 


+ H 2 2 


-> 2 H + OH 


&65 




1.40 x 


10" 


- 12 exp(-2000/T) 




M97 


66 


+ C0 2 - 


-> CO + 2 


^66 




2.80 x 


10" 


- n exp(-26500/T) 




M97 
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Table 1 — Continued 



Number 


Reaction 




Rate Coefficient a 




Reference b 


67 


H+ + -» 


0+ +H 


h 7 = 


7.00 x 10- 10 cxp(-232/T) 




M97 


68 


H 2 + C -» CH 2 + 7 


hs = 


1.00 x 10~ 17 




M97 


69 


H 2 +CHh 


• CH 2 + H 


hg = 


2.38 x 10~ 10 exp(-1760/T) 




M97 


70 


H 2 +CHh 


• C + H + H 2 


ho = 


6.00 x 10~ 9 exp(-40200/T) 




M97 


71 


H 2 + CH 2 - 


-» CH 3 + H 


k 71 = 


5.18 x 10- n (T/300) 17 cxp(- 


•6400/T) 


M97 


72 


H 2 + CH3 - 


-> CH 4 + H 


k 72 = 


6.86 x 10- 14 (T/300) 2 74 cxp(- 


-4740/T) 


M97 


73 


H 2 + OH - 


■ H 2 + H 


hz = 


1.55 x 10- 12 (T/300) 1 - 60 exp(- 


-1660/T) 


M97 


74 


H 2 + OH - 


• O + H + H 2 


k 7i = 


6.00 x 10- 9 exp(-50900/T) 




M97 


75 


H 2 +H 2 0- 


-» OH + H + H 2 


k 7 5 = 


5.80 x 10- 9 cxp(-52900/T) 




M97 


76 


H 2 + 2 -> 


2 H + H 


he = 


2.40 x 10- 10 exp(-28500/T) 




M97 


77 


H 2 + 2 -» 


20H 


k 77 = 


3.16 x 10- 10 exp(-21890/T) 




M97 


78 


H 2 + 2 -» 


20 + H 2 


hs = 


6.00 x 10- 9 cxp(-52300/T) 




M97 


79 


H 2 + 2 H - 


^H 2 2 + H 


k 7 g = 


1.20 x 10- 12 exp(-9400/T) 




M97 


80 


H 2 + C0 2 - 


-> H 2 + CO 


kso = 


2.94 x 10- 14 (T/300)°- 50 exp(- 


-7550/T) 


M97 


81 


H+ + 2 -> 


2 H+ + H 2 


kgi = 


5.00 x 10- 9 cxp(-150/T) 




M97 


82 


C+ + H 2 -» 


CH+ +H 


k&2 = 


1.00 x 10- 10 exp(-4640/T) 




M97 


83 


CH + CH 4 


— > CH 2 + CH 3 


k S 3 = 


2.28 x 10- 11 (T/300) a70 exp(- 


-3000/T) 


M97 


84 


CH + OH - 


■* HCO + H 


fc 84 = 


1.44 x 10- 11 (T/300)°- 50 exp(- 


-5000/T) 


M97 


85 


CH + HCO 


-» CH 2 + CO 


ksb = 


2.87 x 10- 12 (T/300)°- 70 exp(- 


-500/T) 


M97 


86 


CH + H 2 CO -» CH 2 +HCO 


km = 


9.21 x 10- 12 (T/300)°- 70 exp(^ 


-2000/T) 


M97 


87 


CH + 2 - 


■ HCO + O 


h 7 = 


1.44 x 10- 11 (T/300)°- 50 cxp(- 


-3000/T) 


M97 


88 


CH + 2 H 


-» CH 2 + 2 


^88 = 


2.94 x 10- 13 (T/300)°- 50 cxp(- 


-7550/T) 


M97 


89 


CH + 2 H 


-» HCO + OH 


^89 = 


1.44 x 10- 11 (T/300)°- 50 exp(- 


-3000/T) 


M97 


90 


CH + C0 2 


-» HCO + CO 


^90 = 


2.94 x 10- 13 (T/300)°- 50 cxp(- 


-3000/T) 


M97 


91 


2CH 2 -» CH + CH 3 


hi = 


4.00 x 10- 10 cxp(-5000/T) 




M97 


92 


CH 2 + CH4 


^2CH 3 


h 2 = 


7.13 x 10- 12 cxp(-5050/T) 




M97 


93 


CH 2 + HCO -» CH 3 + CO 


hz = 


3.00 x 10" 11 




M97 


94 


CH 2 + H 2 CO -» CH 3 + HCO 


k 94 = 


3.30 x 10- 13 exp(-3270/T) 




M97 


95 


CHj + H - 


> CH+ + H 2 


hb = 


1.00 x 10- 9 exp(-7080/T) 




M97 


96 


CH 3 + H 2 CO -» CH 4 + HCO 


h& = 


2.94 x 10- 13 (T/300) a50 exp(- 


-3000/T) 


M97 


97 


CH+ + H - 


* CH+ + H 2 


h 7 = 


7.00 x 10- 10 exp(-10560/T) 




M97 


98 


OH + CH 2 


-» CH 3 + O 


hs = 


1.44 x 10- 11 (T/300)°- 50 cxp(- 


-3000/T) 


M97 


99 


OH + CH 2 


-» H 2 + CH 


hg = 


1.44 x 10- n (T/300)°' 50 exp(- 


-3000/T) 


M97 


100 


OH + CH 2 


-» H 2 CO + H 


hoo = 


= 3.00 x 10" 11 




M97 


101 


OH + CH 3 


-» H 2 + CH 2 


hoi = 


= 5.42 x 10- n (T/300) - 70 exp(- 


-1000/T) 


M97 


102 


OH + CH 4 


-» H 2 + CH 3 


kl02 = 


= 8.87 x lO-^^T/SOO^-^exp^- 


-1400/T) 


M97 


103 


20H^ H 2 + 


kl03 = 


= 1.65 x 10- 12 (T/300) 1 14 exp( 


-50/T) 


M97 


104 


OH + CO- 


■+ C0 2 + H 


hoi = 


= 1.00 x 10~ 13 




M97 


105 


OH + H 2 2 


; ^H 2 + 2 H 


kl05 = 


= 1.30 x 10- n cxp(-670/T) 




M97 


106 


H 2 + CH 3 


CH 4 + OH 


kl06 = 


= 1.20 x 10- 14 (T/300) 2 - 90 exp( 


-7480/T) 


M97 


107 


CO + 2 - 


>co 2 + o 


ho 7 = 


= 4.20 x 10- 12 exp(-24000/T) 




M97 
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Number Reaction Rate Coefficient a Reference b 



108 


CO + 2 H -> C0 2 + OH 


fcios = 2.50 x 


10" 


■ 10 exp(-11900/T) 




M97 


109 


2 + CH 2 -» HCO + OH 


fciog = 4.10 x 


10" 


- n exp(-750/T) 




M97 


110 


2 + CH 2 -» H 2 CO + O 


feno = 1.44 x 


10" 


11 (T/300)°- 50 exp(- 


3500/T) 


M97 


111 


2 + CH 3 -» 2 H + CH 2 


fcin = 5.30 x 


10" 


12 exp(-34975/T) 




M97 


112 


2 + CH 3 -» H 2 CO + OH 


fcn2 = 5.50 x 


10" 


■ 13 exp(-4500/T) 

*- \ If 




M97 


113 


2 + CH 4 -» CH 3 + 2 H 


fcii3 = 6.60 x 


10" 


- n exp(-28630/T) 




M97 


114 


2 + HCO^ 2 H + CO 


fcn4 = 8.50 x 


10" 


- n exp(-850/T) 




M97 


115 


2 H + CH 3 -» CH 4 + 2 


fcn 5 = 6.00 x 


10" 


-12 




M97 


116 


2 H + H 2 -» H 2 2 + OH 


fcne = 4.60 x 


10" 


- n exp(-16500/T) 




M97 


117 


2 H + HCO^H 2 CO + 2 


feii7 = 5.00 x 


10" 


-11 




M97 


118 


2 H + H 2 CO -» H 2 2 + HCO 


fcn 8 = 3.30 x 


10" 


- 12 exp(-5870/T) 




M97 


119 


20 2 H -> H 2 2 + 2 


fcng = 3.10 X 


10" 


- 12 exp(-775/T) 




M97 


120 


H + HCO^ CO + H 2 


fci 2 o = 1-50 x 


10" 


-10 




M97 


121 


C + H -» CH + 7 


fcl21 = 1.00 X 


10" 


-17 




M97 


122 


2C^C 2 + 7 


fc 122 = 1.00 x 


10" 


-17 




M97 


123 


C + O^ CO + 7 


fc 123 = 2.10 x 


10" 


-19 




M97 


124 


C + CH -» C 2 + H 


fc i24 = 1.73 x 


10" 


-11 




M97 


125 


C + OH^ CO + H 


kl25 = 1.10 x 


10" 


- 10 (T/300) - 50 




M97 


126 


C + HCO^ CH + CO 


fci 26 = 1.00 X 


10" 


-li 




M97 


127 


C + 2 -» CO + O 


fci 27 = 3.30 X 


10" 


-li 




M97 


128 


CH 3 + HCO -» CH 4 + CO 


fci 28 = 2.00 X 


10" 


-10 




M97 


129 


O + H -» OH + 7 


fci 29 = 9.90 X 


10" 


- 19 (T/300)-°- 38 




M97 


130 


20 -» 2 +7 


fci 30 = 4.90 x 


10" 


-^(T/aoo) 1 - 58 




M97 


131 


O + CH^ CO + H 


fci3i = 6.60 x 


10" 


li 




M97 


132 


+ CH^HCO++e 


fc 132 = 2.00 x 


10" 


- n (T/300) - 44 




M97 


133 


O + CH 2 -y CO + 2H 


fci 33 = 2.00 X 


10" 


-10 




M97 


134 


O + CH 3 -» H 2 CO + H 


fc i34 = 1.40 x 


10" 


-10 




M97 


135 


O + OH -» 2 + H 


fc 135 = 4.34 x 


10" 


11 (T/300)-°- 50 exp( 


-30/T) 


M97 


136 


o + c 2 -» co + c 


&136 = 5.00 X 


10" 


- n (T/300) - 50 




M97 


137 


O + HCO^ C0 2 + H 


fci 37 = 5.00 X 


10" 


-li 




M97 


138 


O + HCO^ OH + CO 


^138 = 5.00 X 


10" 


-li 




M97 


139 


O + 2 H -» OH + 2 


fci 39 = 5.40 X 


10" 


-li 




M97 


140 


OH + HCO^H 2 + CO 


k U Q = 1.70 x 


10" 


-10 




M97 


141 


OH + H 2 CO -» H 2 + HCO 


fcui = 4.98 x 


10" 


n (T/300)exp(-601/T) 


M97 


142 


OH + 2 H^H 2 + 2 


fc i42 = 4.00 x 


10" 


-li 




M97 


143 


2HCO -» H 2 CO + CO 


fci 43 = 5.00 X 


10" 


-li 




M97 


144 


H+ + CH -> CH+ + H 


fc i44 = 1.90 x 


10" 


-9 




M97 


145 


H+ + CH 2 -» CH+ + H 2 


fci45 = 1.40 X 


10" 


-9 




M97 


146 


H+ + CH 2 -» CH+ + H 


fci 46 = 1.40 x 


10" 


-9 




M97 


147 


H+ + CH 3 -» CH+ + H 


fci 47 = 3.40 X 


10" 


-9 




M97 


148 


H+ + CH 4 -» CH+ + H 2 


fci 48 = 2.30 X 


10" 


-9 




M97 
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Number Reaction Rate Coefficient a Reference b 



149 


H+ + CH 4 -» CH+ + H 


fcl 4 g = 1.50 X 


10" 


-9 


M97 


150 


H+ + OH — > OH+ + H 


&150 = 2.10 x 


10" 


-9 


M97 


151 


H+ +H 2 -> H 2 0+ + H 


fci5i = 6.90 x 


10" 


-9 


M97 


152 


H+ + C 2 -> C+ + H 


fc 152 = 3.10 x 


10" 


-9 


M97 


153 


H+ + HCO -» CO+ + H 2 


fcl53 = 9.40 X 


10" 


-10 


M97 


154 


H+ + HCO -» H+ + CO 


fci 54 = 9.40 x 


10" 


-10 


M97 


155 


H+ + HCO -» HCO+ + H 


fcl55 = 9.40 X 


10" 


-10 


M97 


156 


H++H 2 CO^H 2 CO+ + H 


fci 56 = 2.96 x 


10" 


-9 


M97 


157 


H+ +H 2 CO -> HCO+ +H 2 


fci57 = 3.57 x 


10" 


-9 


M97 


158 


H+ + 2 -> Oj + H 


fci 58 = 2.00 X 


10" 


-9 


M97 


159 


H+ + C0 2 -» HCO+ + 


ki 59 = 3.50 x 


10- 


-9 


M97 


160 


H- + C^ CH + e 


fci 60 = 1.00 X 


10" 


-9 


M97 


161 


H" + 0^ OH + c 


fciei = 1.00 x 


10" 


-9 


M97 


162 


H-+CH^CH 2 +c 


fcl62 = 1-00 X 


10" 


-10 


M97 


163 


H" + CH 2 -» CH 3 + c 


fcl 63 = 1.00 X 


10" 


-9 


M97 


164 


H- + CH 3 -» CH 4 + c 


fc 164 = 1.00 x 


10- 


-9 


M97 


165 


H" +OH^ H 2 + c 


fci 65 = 1.00 X 


10" 


-10 


M97 


166 


H-+CO^HCO + c 


&166 = 5.00 X 


10" 


-11 


M97 


167 


H-+HCO^H 2 CO + e 


fcl 67 = 1.00 X 


10" 


-9 


M97 


168 


H+ + C -» CH+ + H 


fci 68 = 2.40 X 


10" 


-9 


M97 


169 


H+ + O -» OH+ + H 


fci 6 g = 1.50 X 


10" 


-9 


M97 


170 


H+ + CH -» CH+ + H 2 


fci 70 = 7.10 x 


10" 


-10 


M97 


171 


H+ + CH -> CH+ + H 

z z 1 


k 171 = 7.10 x 


10" 


-10 


M97 


172 


H+ + CH 2 -» CH+ + H 

Z 1 ^ o 


fci 72 = 1.00 x 


10" 


-9 


M97 


173 


+ CH2 — > + H2 


fci 73 = 1.00 X 


10" 


-9 


M97 


174 


H+ + CH 4 -» CH+ + H 2 


fc i74 = 1.40 x 


10" 


-9 


M97 


175 


H+ + CH 4 -> CH5+ + H 


/C175 = 1.14 X 


10" 


-10 


M97 


176 


H+ + CH 4 -» CH+ + H + H 2 


fci 76 = 2.30 X 


10" 


-9 


M97 


177 


H+ + OH -» OH+ + H 2 


fci77 = 7.60 X 


10" 


-10 


M97 


178 


H+ +OH^ H 2 0+ +H 


fcl78 = 7.60 X 


10" 


-10 


M97 


179 


H++H 2 O^H 2 0+ + H 2 


fcl79 = 3.90 X 


10" 


-9 


M97 


180 


H+ + H 2 -» H3O+ + H 


fci 80 = 3.40 x 


10" 


-9 


M97 


181 


H+ + C 2 -> C+ + H 2 


fcisi = 1.10 x 


10" 


-9 


M97 


182 


H++CO^CO++H 2 


fci 82 = 6.40 X 


10" 


-10 


M97 


183 


H++CO^HCO++H 


fc 183 = 2.16 x 


10" 


-9 


M97 


184 


H+ + HCO -» H+ + CO 


fci 84 = 1.00 x 


10" 


-9 


M97 


185 


H+ + HCO -» HCO+ + H 2 


fci 85 = LOO X 


10" 


-9 


M97 


186 


H+ + H 2 CO -» H 2 CO+ + H 2 


fcl86 = 1-40 X 


10" 


-9 


M97 


187 


H+ + H 2 CO -» HCO+ + H + H 2 


fcl87 = 1.40 X 


10" 


-9 


M97 


188 


H+ + 2 -» 2 H+ + H 


kiss = 1-90 x 


10" 


-9 


M97 


189 


H+ + 2 -» 0+ + H 2 


fci 89 = 8.00 X 


10" 


-10 


M97 
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Reaction 


Rate Coefficient 


a 


Reference b 


190 


H+ + C0 2 - 


♦ HCOj + H 


k wo = 2.35 x 


10" 


-9 


M97 


191 


H+ + C -» CH+ + H 2 


fcigi = 2.00 X 


10" 


-9 


M97 


192 


H+ + OH+ + H 2 


fci 92 = 8.00 x 


10" 


-10 


M97 


193 


H+ + CH -> 


CH+ + H 2 


fci 93 = 1.20 X 


10" 


-9 


M97 


194 


H+ + CH 2 - 


» CH+ + H 2 


fc 194 = 1.70 x 


10" 


-9 


M97 


195 


H+ + CH 3 - 


♦ CH+ + H 2 


fc 195 = 2.10 x 


10" 


-9 


M97 


196 


H+ + CH 4 - 


> CH5+ + H 2 


fci 96 = 2.40 X 


10" 


-9 


M97 


197 


Hg + OH — > 


H 2 0+ + H 2 


fci 97 = 1.30 X 


10" 


-9 


M97 


198 


H+ +H 2 - 


* H 3 0+ + H 2 


fc 198 = 5.90 x 


10" 


-9 


M97 


199 


Hg + CO — > 


HCO+ + H 2 


fcigg = 1.73 X 


10" 


-9 


M97 


200 


+ HCO - 


-» H 2 CO+ + H 2 


fc 2 oo = 1-70 x 


10" 


-9 


M97 


201 


H+ +H 2 CO 


-» H 3 CO+ + H 2 


fc 20 i = 6.30 x 


10" 


-9 


M97 


202 


H+ + C0 2 - 


* HCOj + H 2 


fc 202 = 2.00 x 


10" 


-9 


M97 


203 


He+ + CH - 


> C+ + H + He 


fc 20 3 = 1.10 X 


10" 


-9 


M97 


204 


Hc+ + CH - 


> CH+ + He 


fc 2 04 = 5.00 x 


10" 


-10 


M97 


205 


He+ + CH 2 - 


-> CH+ + H + He 


k 2 05 = 7.50 x 


10" 


-10 


M97 


206 


He+ + CH 2 - 


-» C+ + H 2 + He 


&206 = 7.50 X 


10" 


-10 


M97 


207 


He+ + CH 3 - 


-> CH+ + H 2 + He 


k 207 = 1.80 x 


10" 


-9 


M97 


208 


He+ + CH 4 - 


-» CH 3 + H+ + He 


fc 2 08 = 4.80 x 


10" 


-10 


M97 


209 


He+ + CH 4 - 


-» CH+ + H 2 + He + H 


fc 20 g = 2.40 x 


10" 


-10 


M97 


210 


He+ + CH 4 - 


-» CH+ + H 2 + He 


fc 210 = 9.50 x 


10" 


-10 


M97 


211 


Hc+ + CH 4 - 


-» CH+ + He + H 


fc 2 n = 8.50 x 


10" 


-11 


M97 


212 


He+ + CH 4 - 


CH+ + He 


fc 212 = 5.10 x 


10" 


-11 


M97 


213 


He+ + OH - 


* 0+ + H + He 


&213 = 1-10 X 


10" 


-9 


M97 


214 


He+ + H 2 


-» H+ + OH + He 


fc 2i4 = 2.04 x 


10" 


-10 


M97 


215 


He+ + H 2 


-» OH+ + H + He 


k 2 i5 = 2.86 x 


10" 


-10 


M97 


216 


He+ + H 2 


-» H 2 0+ + He 


fei6 = 6.05 x 


10" 


-11 


M97 


217 


He+ + C 2 -» 


C+ + C + He 


fc 2 i 7 = 1.60 x 


10" 


-9 


M97 


218 


Hc+ + C 2 -» 


Cj +Hc 


fc 218 = 5.00 x 


10" 


-10 


M97 


219 


He+ + CO - 


♦ C+ + O + He 


fc 2 ig = 1.60 x 


10" 


-9 


M97 


220 


He+ + HCO 


-» CO+ + H + He 


fc 220 = 4.90 x 


10" 


-10 


M97 


221 


He+ + HCO 


-» CH+ + O + He 


fc 22 i = 4.90 x 


10" 


-10 


M97 


222 


He+ + HCO 


-» HeH+ + CO 


fc 222 = 3.00 x 


10" 


-10 


M97 


223 


He+ + H 2 CO -> CO+ + H 2 + He 


fc 223 = 1.88 x 


10" 


-9 


M97 


224 


He+ + H 2 CO -» HCO+ + H + He 


fc 224 = 1.14 x 


10" 


-9 


M97 


225 


He+ + 2 -» 


0+ + O + He 


fc 225 = 1.00 x 


10" 


-9 


M97 


226 


He+ + 2 -» 


Oj + Hc 


fc 226 = 3.30 x 


10" 


-11 


M97 


227 


He+ + C0 2 


-» 0+ + C + He 


fc 227 = 1.10 x 


10" 


-11 


M97 


228 


Hc+ + C0 2 - 


-» 0+ + CO + He 


fc 228 = 1.00 x 


10" 


-10 


M97 


229 


He+ + C0 2 - 


-» CO+ + O + He 


fc 229 = 8.70 x 


10" 


-10 


M97 


230 


He+ + C0 2 - 


-» C+ + 2 + He 


fc 23 o = 4.00 x 


10" 


-11 


M97 
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231 


C+ + H -> CH+ + 7 


&231 





1. 


.70 


X 


10" 


-17 




M97 


232 


C++O^CO++7 


^232 





2 


.50 


X 


10- 


■18 




M97 


233 


C+ + H- -y H + C 


fe33 





2 


.30 


X 


10- 


7 (T/300)" 


0.50 


M97 


234 


C+ + H 2 -> CH+ + 7 


^234 


— 


4 


.00 


X 


10- 


16 (T/300) 


-0.20 


M97 


235 


C+ + CH -» C+ + H 


^235 





3 


.80 


X 


10- 


■10 




M97 


236 


C+ + CH -» CH+ + C 


&236 


= 


3 


.80 


X 


10- 


■ 10 




M97 


237 


C+ + CH 2 -» CH+ + C 


^237 


— 


5 


,20 


X 


10- 


■10 




M97 


238 


C+ + OH -> CO+ + H 


fc 2 38 





7. 


.70 


X 


10- 


■ 10 




M97 


239 


C+ +H 2 -» HCO+ +H 


fc 2 39 





2 


,70 


X 


10- 


9 




M97 


240 


C+ + HCO -» HCO+ + C 


fc 2 40 





4 


.80 


X 


10- 


■ 10 




M97 


241 


C+ + HCO -» CH+ + CO 


&241 





4 


,80 


X 


10" 


■10 




M97 


242 


C+ +H 2 CO -» CH+ + CO 


^242 





2 


,34 


X 


10- 


9 




M97 


243 


C+ + H 2 CO -» HCO+ + CH 


fc 2 43 





7 


.80 


X 


10- 


■ 10 




M97 


244 


C+ +H 2 CO -> H 2 CO+ + C 


&244 





7. 


.80 


X 


10- 


■10 




M97 


245 


C+ + 2 -» CO+ + 


fc 2 45 





3 


.80 


X 


10- 


■10 




M97 


246 


C+ + 2 -» 0+ + CO 


fc 2 4@ 





6 


.20 


X 


10" 


■ 10 




M97 


247 


C+ + C0 2 -» CO+ + CO 


&247 





1 


.10 


X 


10- 


9 




M97 


248 


CH+ + H -» C+ + H 2 


&248 





7. 


.50 


X 


10- 


10 (T/300) 


-0.50 


M97 


249 


CH+ + C -» Cj + H 


fe 2 49 





1 


.20 


X 


10- 


9 




M97 


250 


CH+ + O -» CO+ + H 


&250 





3 


.50 


X 


10- 


■ 10 




M97 


251 


CH+ + H 2 -» CH+ + H 


&251 





1 


.20 


X 


10- 


9 




M97 


252 


CH+ + CH -» C+ + H 2 


&252 





7. 


,40 


X 


10- 


■10 




M97 


253 


CH+ + OH ^ CO+ + H 2 


fe53 





7. 


.50 


X 


10- 


■ 10 




M97 


254 


CH+ + H 2 -» HCO+ + H 2 


fc 2 54 





2 


.90 


X 


10- 


9 




M97 


255 


CH+ + H 2 -» H 3 0+ + C 


fc 2 55 


= 


5 


.80 


X 


10- 


■ 10 




M97 


256 


CH+ +H 2 -» H 2 CO+ +H 


&256 





5 


.80 


X 


10- 


■ 10 




M97 


257 


CH+ + CO -» HCO+ + C 


fe57 


— 


7. 


.00 


X 


10- 


■ 12 




M97 


258 


CH+ + HCO -» CH+ + CO 


&258 





4 


.60 


X 


10- 


-10 




M97 


259 


CH+ + HCO -» HCO+ + CH 


&259 


— 


4 


,60 


X 


10- 


■10 




M97 


260 


CH++H 2 CO^CH+ + CO 


&260 




9 


.60 


X 


10- 


■10 




M97 


261 


CH+ + H 2 CO -» HCO+ + CH 2 


^261 




9 


.60 


X 


10- 


■ 10 




M97 


262 


CH+ + H 2 CO -» H 3 CO+ + C 


&262 




9 


.60 


X 


10- 


■ 10 




M97 


263 


CH+ + 2 -» HCO+ + O 


&263 




9 


.70 


X 


10- 


■ 10 




M97 


264 


CH+ + 2 -» HCO + 0+ 


&264 




1 


.00 


X 


10- 


■ 11 




M97 


265 


CH+ + 2 -» CO+ + OH 


&265 




1 


.00 


X 


10- 


■ 11 




M97 


266 


CH+ + C0 2 HCO+ + CO 


&266 




1 


.60 


X 


10- 


9 




M97 


267 


CH+ + O -» HCO+ + H 


&267 




7. 


.50 


X 


10- 


■10 




M97 


268 


CH+ + H 2 -» CH+ + H 


^268 




1 


.60 


X 


10- 


9 




M97 


269 


CH+ + H 2 -> H3CO+ + H 


&269 




1 


.20 


X 


10- 


9 




M97 


270 


CH+ + HCO -» CH+ + CO 


&270 




4 


.50 


X 


10- 


■10 




M97 


271 


CH+ + H 2 CO -» HCO+ + CH 3 


&271 




2 


.81 


X 


10- 


9 




M97 
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Reaction 



Rate Coefficient a 



Reference b 



272 
273 
274 
275 
276 
277 
278 
279 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
291 
292 
293 
294 
295 
296 
297 
298 
299 
300 
301 
302 
303 
304 
305 
306 
307 
308 
309 
310 
311 
312 



CH+ + 2 



CH+ 
CH+ 
CH+ 
CH+ + H 2 
CH+ 
CH+ 



► HCO+ + OH 
-» H 2 CO+ + CO 
HCO+ + H 2 



C0 2 
O -. 

O^H 2 CO+ + H 



CH5+ + 7 
+ H 2 CO+ + H 2 
■ CH+ + CO 
• HCO+ + CH 3 



OH- 
HCO 
CH+ + HCO 
CH+ + H 2 CO -» HCO+ + CH 4 



CH+ 



H 3 CO+ + 
H+ + 
H + O 
> OH+ + 
-» CH+ 4 
^CO+H 



2 

0++H^ 
0+ +H- - 
0+ + H 2 - 
0+ + CH - 
0+ + CH - 
0+ + CH 2 
0+ + CH 4 
0+ + CH 4 
0+ + OH 
0+ + OH 
0+ +H 2 

0+ + c 2 - 

0+ + C 2 - 
0+ + HCO - 
0+ + HCO 
0+ + H 2 CO 
0+ +H 2 CO 
0+ + 2 
0+ + C0 2 - 
CH+ + H -> 
CH+ + O - 
CH+ + H 2 - 
CH+ + CH 4 
CH++H 2 
CH+ + CO - 

CH++H 2 CO^H 2 CO+ 



H 

O 
H 

■+ CH+ + O 
■* CH+ + OH 
-> CH+ + O 
• OH+ + O 
-0++H 
-» H 2 0+ + O 
C+ + O 
CO+ + c 
-» OH+ + CO 
-» HCO+ + O 
■* HCO+ + OH 
■* H 2 CO+ + O 
O+ + O 



0+ 



CO 

- H9 



CH+ 

CH+ + OH 
> CH5+ + H 
-» CH5+ + CH 3 
-> H3O+ + CH 3 



CH++H 2 CO 



HCO+ + CH 3 
H 2 CO+ 
H 3 CO+ 



CH 4 



-4 

CH 3 



CH+ + 2 -» 0+ + CH 4 



CH+ 



4 T yJ2 —> w 2 

. 4 -r C0 2 -» HCO+ + CH 3 
OH+ + C -> CH+ + O 
OH+ + O -» 0+ + H 



&272 
&273 
&274 
&275 
&276 
fc 2 77 
&278 
fc 2 79 
&280 
&281 
&282 
^283 
&284 
&285 
^286 
&287 
&288 
&289 
&290 
&291 
&292 
^293 
&294 
&295 
^296 
&297 
&298 
^299 
&300 
&301 
&302 
^303 
&304 
&305 
&306 
&307 
&308 
&309 
&310 
fell 
fel2 



9.10 x 10 
: 1.60 x 10 
: 4.00 x 10 
: 4.00 x 10 
: 1.30 x 10 
7.20 x 10 
4.40 x 10 
4.40 x 10 
: 1.60 x 10 
5.00 x 10 
6.80 x 10 
: 2.30 x 10 
1.70 x 10 
3.50 x 10 
3.50 x 10 
9.70 x 10 
1.10 x 10 
8.90 x 10 
3.60 x 10 
3.60 x 10 
3.20 x 10 
4.80 x 10 
4.80 x 10 
4.30 x 10 
4.30 x 10 
1.40 x 10 
2.10 x 10 
1.90 x 10 
: 9.40 x 10 
: 1.00 x 10 
: 1.00 x 10 
3.30 x 10 
1.50 x 10 
2.60 x 10 
: 1.40 x 10 
: 1.62 x 10 
1.98 x 10 
4.00 x 10 
: 1.20 x 10 
1.20 x 10 
7.10 x 10 



-10 
-9 
-10 
-11 

-14 

-10 

-10 

-10 

-9 

-12 

-10 



(T/300)- 1 



" 7 (T/300) 



-0.50 



-9 
•10 
-10 
-10 
-10 
-10 
-10 
-10 
9 

-10 
-10 
-10 
-10 
-9 
-9 
-11 
-10 
-11 
9 

-11 

-9 

-9 

-9 

-9 

-9 

-10 



-10 



M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
M97 
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Reaction 


Rate CoefBcicnt 


a 


Reference b 


313 


OH+ + H 2 -> 


H 2 0+ + H 


fc 3 i 3 = 1.01 x 


10" 


-9 


M97 


314 


OH+ + CHh 


CH+ + OH 


k 31 4 = 3.50 x 


10" 


-10 


M97 


315 


OH+ + CHh 


CH+ + 


k 315 = 3.50 x 


10" 


-10 


M97 


316 


OH+ + CH 2 - 


-> CH+ + OH 


fc 316 = 4.80 x 


10" 


-10 


M97 


317 


OH+ + CH 2 - 


■* CH+ + 


fc 317 = 4.80 x 


10" 


-10 


M97 


318 


OH+ + CH 4 - 


-» H 3 0+ + CH 2 


&318 = 1-31 X 


10" 


-9 


M97 


319 


OH+ + CH 4 - 


-» CH5+ + 


fc3iq = 1.95 X 


10" 


-10 


M97 


320 


OH+ + OH - 


■ H 2 0+ + 


^320 = 7.00 X 


10" 


-10 


M97 


321 


OH+ +H 2 - 


-» H 2 0+ + OH 


A: 321 = 1.59 x 


10" 


-9 


M97 


322 


OH+ + H 2 - 


-» H 3 0+ + O 


fc322 = 1-30 X 


10" 


-9 


M97 


323 


OH+ + C 2 -> 


Cj + OH 


fc323 = 4.80 X 


10" 


-10 


M97 


324 


OH+ + CO -> 


• HCO+ + O 


&324 = 1-05 X 


10" 


-9 


M97 


325 


OH+ + HCO 


-» H 2 0+ + CO 


k 3 25 = 2.80 X 


10" 


-10 


M97 


326 


OH+ + HCO 


-» HCO+ + OH 


fc 326 = 2.80 x 


10" 


-10 


M97 


327 


OH+ + HCO 


-> H 2 CO+ + 


fc 327 = 2.80 x 


10" 


-10 


M97 


328 


OH+ + H 2 CO 


-» H3CO+ + O 


&328 = 1-12 X 


10" 


-9 


M97 


329 


OH+ + H 2 CO 


-» H 2 CO+ + OH 


^329 = 7.44 X 


10" 


-10 


M97 


330 


OH+ + 2 -> 


Oj + OH 


fc 330 = 5.90 x 


10" 


-10 


M97 


331 


OH+ + C0 2 - 


-» HCO+ + O 


/c 331 = 1.44 x 


10" 


-9 


M97 


332 


CH5+ + H -> 


CH+ + H 2 


fc 332 = 2.00 x 


10" 


-11 


M97 


333 


CH5+ + C -» 


CH+ + CH 4 


fc 333 = 1.20 x 


10" 


-9 


M97 


334 


CH5+ + -» 


H3O+ + CH 2 


fc 334 = 2.20 x 


10" 


-10 


M97 


335 


CH5+ + -» 


H 3 CO+ + H 2 


fc 335 = 4.40 x 


10" 


-12 


M97 


336 


CH5+ + CH - 


CH+ + CH 4 


fc 336 = 6.90 x 


10" 


-10 


M97 


337 


CH5+ + CH 2 


— > CHg^ + CH 4 


fc 337 = 9.60 x 


10" 


-10 


M97 


338 


CH5+ + OH - 


-» H 2 0+ + CH 4 


^338 = 7.00 X 


10" 


-10 


M97 


339 


CH5+ +H 2 


-» H3O+ + CH 4 


^339 = 3.70 X 


10" 


-9 


M97 


340 


CH5+ + CO - 


-> HCO+ + CH 4 


fc 340 = 1.00 x 


10" 


-9 


M97 


341 


CH5+ + HCO 


-» H 2 CO+ + CH 4 


/c 341 = 8.50 x 


10" 


-10 


M97 


342 


CH5+ + H 2 CO -» H 3 CO+ + CH 4 


&342 = 4.50 X 


10" 


-9 


M97 


343 


CH5+ + C0 2 


-» HCO+ + CH 4 


fc 343 = 3.20 x 


10" 


-11 


M97 


344 


H 2 0+ + C^ 


CH+ + OH 


&344 = 1-10 X 


10" 


-9 


M97 


345 


H 2 0+ + -> 


0++H 2 


fc 345 = 4.00 x 


10" 


-11 


M97 


346 


H 2 0++H 2 - 


>H 3 0++H 


fc 346 = 8.30 x 


10" 


-10 


M97 


347 


H 2 0+ + CH- 


■* CH+ + H 2 


fc 347 = 3.40 x 


10" 


-10 


M97 


348 


H 2 0+ + CH- 


■* CH+ + OH 


fc 348 = 3.40 x 


10" 


-10 


M97 


349 


H 2 0+ + CH 2 


-> CH+ + OH 


fc 349 = 4.70 x 


10" 


-10 


M97 


350 


H 2 0+ + CH 2 


-» CH+ +H 2 


fcsso = 4.70 x 


10" 


-10 


M97 


351 


H 2 0+ + CH 4 


-» H3O+ + CH 3 


fc 35 i = 1.11 X 


10" 


-9 


M97 


352 


H 2 0+ + OH - 


-+ H3O+ + O 


&352 = 6.90 X 


10" 


-10 


M97 


353 


H 2 0++H 2 


-» H3O+ + OH 


/C353 = 2.10 x 


10" 


-9 


M97 
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Table 1 — Continued 



Number 


Reaction 






Rate Coefficient a 


Reference 


354 


H 2 0+ + C 2 ^ C^ + H 2 


^354 


= 


4.70 x 


10" 


-10 


M97 


355 


H 2 0+ + CO -» HCO+ + OH 


^355 


= 


5.00 x 


10" 


-10 


M97 


356 


H 2 0+ + HCO -» H 3 0+ + CO 


^356 


= 


2.80 x 


10" 


-10 


M97 


357 


H 2 0+ + HCO -» HCO+ + H 2 


^357 


= 


2.80 x 


10" 


-10 


M97 


358 


H 2 0+ + HCO -» H 2 CO+ + OH 


&358 


= 


2.80 x 


10" 


-10 


M97 


359 


H 2 0+ + H 2 CO -» H 2 CO+ + H 2 


^359 


= 


1.41 x 


10" 


-9 


M97 


360 


H 2 0+ + H 2 CO -» H 3 CO+ + OH 


^360 


= 


6.62 x 


10" 


-10 


M97 


361 


H 2 0+ + 2 -» 0+ +H 2 


^361 


= 


4.30 x 


10" 


-10 


M97 


362 


H3O+ + C -» HCO+ + H 2 


h62 


= 


1.00 x 


10" 


-11 


M97 


363 


H 3 0+ + H" -» OH + H 2 + H 


h63 


= 


2.30 x 


10" 


- 7 (T/300)-°- 50 


M97 


364 


H 3 0+ + H- -» H 2 + H 2 


&364 


= 


2.30 x 


10" 


- 7 (T/300)-°- 50 


M97 


365 


H 3 0+ + CH^CH+ + H 2 


&365 


= 


6.80 x 


10" 


-10 


M97 


366 


H 3 0+ + CH 2 -» CH+ + H 2 


^366 


= 


9.40 x 


10" 


-10 


M97 


367 


H 3 0+ + H 2 CO -» H 3 CO+ + H 2 


&367 


= 


3.40 x 


10" 


-9 


M97 


368 


c+ + c -» c+ + c 2 


^368 


= 


1.10 x 


10" 


-10 


M97 


369 


C+ + O -» CO+ + c 


&369 


= 


3.10 x 


10" 


-10 


M97 


370 


C+ + CH -» CH+ + C 2 


&370 


= 


3.20 x 


10" 


-10 


M97 


371 


C^ + CH 2 — ► CH^ + C 2 


&371 


= 


4.50 x 


10" 


-10 


M97 


372 


C+ + OH -» OH+ + C 2 


&372 


= 


6.50 x 


10" 


-10 


M97 


373 


C+ +HCO -» HCO+ + C 2 


&373 


= 


3.80 x 


10" 


-10 


M97 


374 


c+ + 2 -» CO+ + CO 


&374 


= 


8.00 x 


10" 


-10 


M97 


375 


CO++H^H+ + CO 


&375 


= 


7.50 x 


10" 


-10 


M97 


376 


CO+ + C^C+ + CO 


&376 


= 


1.10 x 


10" 


-10 


M97 


377 


CO+ + 0^0+ + CO 


&377 


= 


1.40 x 


10" 


-10 


M97 


378 


CO+ +H 2 -» HCO+ +H 


&378 


= 


1.80 x 


10" 


-9 


M97 


379 


CO+ + CH -» CH+ + CO 


&379 


= 


3.20 x 


10" 


-10 


M97 


380 


CO+ + CH -» HCO+ + C 


^380 


= 


3.20 x 


10" 


-10 


M97 


381 


CO+ + CH 2 -» CH+ + CO 


^381 


= 


4.30 x 


10" 


-10 


M97 


382 


CO+ + CH 2 -» HCO+ + CH 


&382 


= 


4.30 x 


10" 


-10 


M97 


383 


CO+ + CH 4 -» CH+ + CO 


&383 


= 


7.93 x 


10" 


-10 


M97 


384 


CO+ + CH 4 -» HCO+ + CH 3 


&384 




4.55 x 


10" 


-10 


M97 


385 


CO+ + OH -» OH+ + CO 


&385 




3.10 x 


10" 


-10 


M97 


386 


CO+ + OH -» HCO+ + O 


^386 




3.10 x 


10" 


-10 


M97 


387 


CO+ +H 2 -» H 2 0+ + CO 


hs7 




1.72 x 


10" 


-9 


M97 


388 


CO+ +H 2 -» HCO+ + OH 


hs8 




8.84 x 


10" 


-10 


M97 


389 


CO+ + C 2 -» C+ + CO 


^389 




8.40 x 


10" 


-10 


M97 


390 


CO+ + HCO -» HCO+ + CO 


^390 




7.40 x 


10" 


-10 


M97 


391 


CO+ + H 2 CO -» HCO+ + HCO 


^391 




1.65 x 


10" 


-9 


M97 


392 


CO+ + H 2 CO -» H 2 CO+ + CO 


&392 




1.35 x 


10" 


-9 


M97 


393 


CO+ + 2 -» 0+ + CO 


&393 




1.20 x 


10" 


-10 


M97 


394 


HCO+ + C -» CH+ + CO 


&394 




1.10 x 


10" 


-9 


M97 



- 23 - 



Table 1 — Continued 



Number 


Reaction 






Rate Coefficient a 


Rclerence 


395 


HCO+ 4- H — > CO 4- Ho 

ll\Jvy | 11 7 V Vy \ 11^ 


^395 




2.30 x 


10" 


" 7 (T/300)-°- 50 


M97 


396 


HCO+ 4- CH — » CHt 4- CO 


ft*ooc 

^396 




6.30 x 


10" 


-10 


M97 


397 


HCO+ 4- CHo — > OTT+ 4- CO 

1 1 v. V y \ Vyll^ 7 Vyllo \^ Vj\y 


^397 




8.60 x 


10" 


-10 


M97 


398 


HCO+ 4- OH — > HoO+ + CO 

1 1 v vy i vyii 7 2 y^ v_ v ./ 


i- 

"-398 




6.20 x 


10" 


-10 


M97 


399 


HCO+ 4- OH — » HCO+ 4- H 

llV^Vy | Vyll ' 1 1 V. V / i | -L.L 


^399 




1.00 x 


10" 


-9 


M97 


400 


HCO+ 4- HoO — > HqO+ 4- CO 

x x v_ / vy | iiy^ vy ' iij vy | v_/ vy 


'MOO 




2.50 x 


10" 


-9 


M97 


401 


HCO+ 4- HCO — » H?CO+ + CO 


^40 1 




7.30 x 


10" 


-10 


M97 


402 


HCO+ 4- HtCO — > HqCO+ 4- CO 

1 1 v vy | i iv v.. v / ' i i yj v_ vy \ v_/ vy 


^402 




3.30 x 


10" 


-9 


M97 


403 


HnCO+ 4- CH — » CH+ 4- HoCO 

ll^VyVy | Vjll 7 V 1 1 | 1 ' yj V V7 


^403 




3.10 x 


10" 


-10 


M97 


404 


HoCO+ 4- CH — > CHo" 4- HCO 


A '404 




3.10 x 


10" 


-10 


M97 


405 


HoCO+ 4- CHo — > CHt 4 HCO 

ll^J V.-' Vy | VJ112 7 VyllQ |^ 1 1 V. V y 


^405 




4.30 x 


10" 


-10 


M97 


406 


HoCO+ 4- CHo — > CH+ 4 HoCO 

112 Vy V y | Vyll^ 7 V_yllo |^ ll^VyVy 


^406 




4.30 x 


10" 


-10 


M97 


407 


HoCO+ 4- CH/ — > HqCO+ 4- CHq 

a ±2 -- vy | vyiizj. 7 o ^- vy ^— ^-'--^o 


^407 




9.35 x 


10" 


-11 


M97 


408 


HoCO+ 4- HoO — » HqO+ + HCO 


^408 




2.60 x 


10" 


-9 


M97 


409 


HnCO+ 4- HCO — ► HCO+ + HnCO 

i.i. t 2 i \j\j \ 1 1 v. v / 7 i i v V7 x x 2 v_- vy 


^409 




3.60 x 


10" 


-10 


M97 


410 


HoCO+ 4- HCO — > HqCO+ 4- CO 


Km 

"■410 




3.60 x 


10" 


-10 


M97 


411 


HoCO+ 4- HoCO — > Hc>CO + + HCO 

ii^vVy | l 17 v_ Vy ' l l y v_/ Yy | iiw 






3.20 x 


10" 


-9 


M97 


412 


HoCO+ 4- Oo — > HCO+ 4- OoH 

iij vvy \ vy 2, i i v. \y i vyy,- 1 -- 1 - 


^412 




7.70 x 


10" 


XI 


M97 


413 


H^CO+ 4- CH — ► 0H„" 4- HoCO 

11,} v.. V. y | Vjll 7 Vllo | 1 ' y. V Vy 


^413 




6.20 x 


10" 


-10 


M97 


414 


HqCO+ 4- HoO — > HqO+ + HoOO 

11,} v.. v y lljv ~ lljvy i iijvyvy 7 


^414 




2.30 x 


10" 


-10 


MQ7 


415 


0+ 4- C — > CO+ 4- o 

VJo | V- ~ v Vy | V y 


^415 




5.20 x 


10" 


-11 


M97 


416 


n+ 4- C — » C+ 4- Oo 

vy 2 y^ v ' vy vy 2 


^416 




5.20 x 


10" 


XI 


M97 


417 


+ 4- CH — > CH+ 4- Oo 

Vo I Vyll r V^ll | Vy 2 


^417 




3.10 x 


10" 


-10 


M97 


418 


+ 4- CH — » HCO+ 4 O 

V y y V_/ 1 1 " 1 1 V_ V7 |^ V y 


^418 




3.10 x 


10" 


-10 


M97 


419 


n+ 4- CHo — ► HoCO+ 4 O 

vy |2 '^-'-■'-2 2 ^- ^ 1^ vy 


"-419 




4.30 x 


10" 


-10 


M97 


420 


+ 4- CHo — > CH + 4 Oo 

vy 2 ~ 7 vy±±2 1 vy 2 


^420 




4.30 x 


10" 


-10 


M97 


421 


+ 4- Co — » CO+ 4- CO 

V y y 1 Vy^ 7 v.. v.y VyVy 


^421 




4.10 x 


10" 


-10 


M97 


422 


+ 4- Co — > C + 4 Oo 


^422 




4.10 x 


10" 


-10 


M97 


423 


+ 4- HCO — > OoH+ 4- OO 

Wo 1 1 v vy * vyj AA y^ v_ ■ vy 


&423 




3.60 x 


10" 


-10 


MQ7 


424 


+ 4- HOO — ► HOO+ + Oo 

V y y 1 1 V. V / 7 1 1 V. V / \ Vy 2 


^424 




3.60 x 


10" 


-10 


MQ7 


425 


0+ 4- HoCO — ► HCO+ + Oo 4- H 

V y y 1 1 y V_/ v.y ? 1 1 V_ V y V y y 11 


^425 




2.30 x 


10" 


-10 


M97 


426 


0++H 2 CO^H 2 CO + 4-02 


^426 




2.07 x 


10" 


-9 


M97 


427 


2 H+ + C -y CH+ + 2 


^427 




1.00 x 


10" 


-9 


M97 


428 


2 H+ + -» OH+ + 2 


^428 




6.20 x 


10" 


-10 


M97 


429 


2 H+ + H 2 -> H+ + 2 


^429 




6.40 x 


10" 


-10 


M97 


430 


2 H+ + CH -» CH+ + 2 


^430 




6.20 x 


10" 


-10 


M97 


431 


2 H+ + CH 2 -» CH+ + 2 


^431 




8.50 x 


10" 


-10 


M97 


432 


2 H+ + OH -> H 2 0+ 4- 2 


^432 




6.10 x 


10" 


-10 


M97 


433 


2 H+ + H 2 -» H 3 0+ + 2 


^433 




8.20 x 


10" 


-10 


M97 


434 


2 H+ + CO -» HCO+ + 2 


^434 




8.40 x 


10" 


-10 


M97 


435 


2 H+ + HCO -» H 2 CO+ 4- 2 


^435 




7.10 x 


10" 


-10 


M97 



Table 1 — Continued 



Number 



Reaction 



Rate Coefficient a 



Reference b 



436 
437 
438 
439 
440 
441 
442 
443 
444 
445 
446 
447 
448 
449 
450 
451 
452 
453 
454 
455 
456 
457 
458 
459 
460 
461 
462 
463 
464 
465 
466 
467 
468 
469 
470 
471 
472 
473 
474 
475 
476 



2 H+ + H 2 CO -» H 3 CO+ + 2 
2 H+ + C0 2 -» HCO+ + 2 
HCO+ + C -» CH+ + C0 2 
HCO+ + O -» HCO+ + 2 
HCO+ + CH 4 -» CH5+ + C0 2 
HCO+ + H 2 -» H 3 0+ + C0 2 
HCO+ + CO -» HCO+ + C0 2 
C++e^C + 7 
CH+ + c -» C + H 
-» C + H 2 
CH + H 
-» CH 2 + H 
■+ CH + H 2 
-> CH + 2H 
- CH 3 + 7 
+ 7 

CH 3 + H 
-> CH 2 + 2H 

O + H 
— > CH 4 + H 
— > CH 3 + H 2 
-» OH + H 
-> + H 2 
^H 2 + H 
-» OH + 2H 
2C 

-» o + c 

CO + H 
HCO + H 
CO + 2H 
H 2 CO + 7 
CO + H + H 2 
HCO + 2H 
H 2 CO + H 



CE£ + c - 
CE.% + c - 
CH+ + e - 
CH+ + c - 
CH+ + c - 
CH+ + e - 
0++e^ 
CH+ + c - 
CH+ + c - 
OH+ + c - 
CH5+ +e 
CH5+ +e 
H 2 0++c 
H 2 0+ + c 
H 3 0+ + c 
H 3 0+ + c 
C++c^ 
CO++c^ 
HCO+ + e - 
H 2 CO++c 
H 2 CO++c 
H 2 CO++c 
H 3 CO+ 
H 3 CO+ 
H 3 CO+ 



0++c 



c — » 
c — > 
c — ► 
20 



2 H+ 
HCO+ 
HCO+ 
OH 
H+ 



3^0 2 +H 

C0 2 + H 
e -» CO + OH 
CH 3 -» CH 4 + O 
H 2 CO -» CO+ + H 2 + H 



He+ + C -» C+ + He 



&436 


= 


9.80 x 


10" 


i n 

- 1U 




M97 


^437 


= 


1.10 x 


10" 


-9 




M97 


&438 


= 


1.00 x 


10" 


Q 

-y 




M97 


^439 


= 


1.00 x 


10" 


-9 




M97 


&440 


= 


7.80 x 


10" 


-10 




M97 


fc 44 l 


= 


2.30 x 


10" 


-9 




M97 


&44 2 


= 


7.80 x 


10" 


-If) 




M97 


&443 


= 


4.40 x 


10" 


_12 (T/300) 


— U.Dl 


M97 


&444 


= 


1.50 x 


10" 


" 7 (T/300) _ 


-0 42 


M97 


&445 


= 


1.25 x 


10" 


" 7 (T/300) _ 


-u. ou 


M97 


&446 


= 


1.25 x 


10" 


" 7 (T/300) _ 


-0 50 


M97 


fc 447 


= 


7.75 x 


10" 


" 8 (T/300) _ 


-u.ou 


M97 


fc 44 8 


= 


1.95 x 


10" 


" 7 (T/300) _ 


50 


M97 


fc 44 g 


= 


7.75 x 


10" 


" 8 (T/300)~ 


-0 50 


M97 


&450 


= 


1.10 x 


10" 


" 10 (T/300) 


— 50 


M97 


&451 


= 


3.40 x 


10" 


12 (T/300) 


— 63 


M97 


&452 


= 


1.75 x 


10" 


" 7 (T/300) _ 


50 


M97 


&453 


= 


1.75 x 


10" 


" 7 (T/300) _ 


-u. ou 


M97 


&454 


= 


3.75 x 


10" 


~ 8 (T/300)~ 


u.ou 


M97 


&455 


= 


5.50 x 


10" 


" 7 (T/300) _ 


u.ou 


M97 


&456 


= 


5.50 x 


10" 


7 (T/300)"" 


-0 30 


M97 


&457 


= 


1.60 x 


10" 


" 7 (T/300)~ 


-0 50 


M97 


&458 


= 


2.00 x 


10" 


_7 (T/300)~ 


-0 50 


M97 


&459 


= 


3.50 x 


10" 


" 7 (T/300) _ 


u.ou 


M97 


&460 


= 


6.50 x 


10" 


" 7 (T/300)~ 


-0 50 


M97 


&461 


= 


3.00 x 


10" 


" 7 (T/300)"~ 


-0 50 


M97 


&462 


= 


1.80 x 


10" 


'(T/300) 


-0.50 


M97 


&463 


= 


1.10 x 


10" 


7 (T/300)- 


- 1 .00 


M97 


&464 




1.00 x 


10" 


" 7 (T/300)~ 


-0 50 


M97 


&465 




5.00 x 


10" 


" 7 (T/300)- 


-0.50 


M97 


^466 




1.10 x 


10" 


10 (T/300) 


-0.70 


M97 


&467 




2.00 x 


10" 


" 7 (T/300)- 


-0.50 


M97 


^468 




2.00 x 


10" 


" 7 (T/300)- 


-0.50 


M97 


&469 




2.00 x 


10" 


" 7 (T/300)" 


-0.50 


M97 


&470 




1.95 x 


10" 


- 7 (T/300)" 


-0.70 


M97 


fc 47 l 




3.00 x 


10" 


" 7 (T/300)" 


-0.50 


M97 


k±72 




2.40 x 


10" 


" 7 (T/300)" 


-0.50 


M97 


&473 




1.00 x 


10" 


" 7 (T/300)" 


-0.50 


M97 


fc 474 




3.22 x 


10" 


14 (T/300) 


2 - 20 exp(-2240/T) 


M97 


&475 




1.06 x 


10" 


-9 




M97 


&476 




6.30 x 


10" 


" 15 (T/300) 


0.75 


M97 
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Fig. 1. — The evolutionary trajectories of protostellar clouds with different mctallicitics. The solid curves 
show the temperature of the clouds with metallicity (from the top to bottom) z = 0, 1CP 6 , 1CP 4 , 10~ 2 , 1. The 
dotted curves show the grain temperature for the same clouds. The initial condition is T = 100K, n — 
lcm -3 , y(e) = 1 x 10 -4 , y(H2) = 1 x 10~ 6 . Initially, all the carbon is assumed to be in the form of CII, while 
oxygen is 01. The dashed line marked "tj = 1" represents the line L in the text (Eq.35). 
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Fig. 2. — Initial condition dependence of the evolutionary trajectories. The solid curves are the same as 
those of Figure 1 ("fiducial model"). The initial temperature of dotted curves is 1000K and of the dashed 
curves is 10K, while other quatities are the same as the fiducial model. For same values of metallicity, a 
wide variation in temperature yields convergent evolutionary tracks. 
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Fig. 3. — The contributions to the cooling/heating rates by major processes for the same clouds as illustrated 
in Figure 1: (a)z = 0, (b)z = 1CT 6 , (c)z = 1CT 4 , (d)z = 1CT 2 , and (e)z = 1. Curves labelled "comp", "grain", 
"H 2 form/diss", and "H 2 cont" show the cooling/heating rate by the adiabatic compression, the gas-grain 
heat transfer, formation/dissociation of H 2 , and H 2 continuum. The other curves labelled "H 2 " , "OH+H 2 0" , 
"CO", "CI", and "01" illustrate the line cooling rate by these atoms/molecules. Irregularities in "grain" 
(in panels b-e) and "H 2 " and "OH+H 2 0"(in panel b and c) cooling rates near n ~ 10 14 cm~ 3 are caused by 
sublimaiton of each grain component and subsequent change of the dust opacity. 
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Fig. 4. — The concentrations of hydrogen molecules (solid line), atoms (dotted line) and ions (dashed line) 
for the same clouds as illustrated in Fig.l. 



